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Abstract. The extragalactic background light (EBL) is of fundamental importance both for un- 
derstanding the entire process of galaxy evolution and for 7-ray astronomy. However, the overall 
spectrum of the EBL between 0.1 and 1000 ^m has never been determined directly neither from 
observed luminosity functions (LFs), over a wide redshift range, nor from any multiwavelength 
observation of galaxy spectral energy distributions (SEDs). The evolving, overall spectrum of 
the EBL is derived here utilizing a novel method based on observations only. It is emphasized 
that the local EBL seems already well constrained from the UV up to the mid-IR. Since different 
independent methodologies such as direct measurement, galaxy counts, 7-ray attenuation and 
realistic EBL modelings point towards the same EBL intensity level. A relevant contribution 
from Pop III stars to the local EBL seems unlikely. 
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1. Introduction 

The extragalactic background light (EBL) is the accumulated radiation in the uni- 
verse from star formation process, plus a contribution from active galactic nuclei (AGNs). 
These photons mostly lie in the range ~0. 1-1000 /im. The direct measurement of the EBL 
is a very difficult task subject to high uncertainties. This is mainly due to the contribution 
of zodiacal light, some orders of magnitude larger than the EBL (e.g., Hauser & Dwek 



2001; Chary & Pope 2010). Interestingly, it has been recently claimed by Matsuoka et al. 



(2011) the detection of the EBL free of zodiacal light. Other observational approaches 



set reliable lower limits on the EBL, such as measuring the integrated light from discrete 
extragalactic sources (e.g., Madau & Pozzetti 2000; Fazio et al. 2004; Keenan et al. 2010). 
On the other hand, there are phenomenological approaches in the literature that predict 
an overall EBL model (i.e., between 0.1 and 1000 /zm and for any redshift). These are 
basically of the four kinds described in Dommguez et al. (2011a) and enumerated here in 
Table [1] Generally, any EBL modeling is built from two main quantities: one describing 
the galaxy density evolution over time and another one describing the overall galaxy 
emission (stellar component plus absorption/re-emission by dust). Table [1] briefly sum- 
marizes how these two main quantities are treated in the most relevant EBL modelings 
among the bibliography. We stress that the previous four-types classification is based 
upon how the different methodologies describe the galaxy density evolution. 



We consider the theoretical approach taken in Somerville et al. (2011) Gilmore et al. 



(2011) as complementary to our observationally motivated one to eventually reach a com- 



plete understanding of galaxy evolution. Approaches type (ii) are potentially problematic 
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Type of modeling & refs. 


Galaxy density evolution 


Galaxy emission 


(ij i'brward evolution, 


Semi-analytical mod- 
els 


Modeled. Stellar emission: 


e.g., ISomerville et al. 


IBruzual & Chariot (2003) | (BC03); 


(2011) 


Gilmore et al. 


Dust absorption: |Charlot & Fall| 


JWTj 




1(2000)1 Dust re-emission: templates 






by Rieke et al. (2009) 


(ii) Backward evolution, 
e.g., iFranccschini et al.| 


Observed local-optical 
galaxy luminosity func- 
tion (LF, for starburst 
population) and near-IR 
galaxy LF observed up 
to z — 1.4 (for elliptical 
and spiral populations) 


Modeled. A few galaxy types mor- 
phologically classified based on op- 
tical images. 


(2008) 




(iii) Inferred evolution, 
e.g.,|Finke et al. (2010) 


Parameterization of 

the history of the star 
formation rate density 
of the universe 


Modeled. Stellar emission: Sin- 
gle bursts of solar metallicity from 


Kneiske & Dole (2010) 


Bruzual & Chariot (1993) (in 




Kneiske & Dole 2010)/BC03 (in 
Finke et al. 2010); Dust absorp- 
tion: General extinction law; Dust 
re-emission: Modified black bodies. 
AGN galaxies are not considered. 


(iv) Observed evolution, 
Dommguez et al. (2011a) 


Observed near-IR 
galaxy LF up to z = 4 


Observed. Based on multiwave- 
length photometry from the UV up 
to MIPS 24 for ~ 6000 galaxies up 
to z=l. Consider 25 different galaxy 
types including AGN galaxies. 



Table 1. Classification and comparison of the main characteristics of recent EBL modelings. 



because they imply extrapolations backwards in time of local or low-redshift luminosity 
functions (LFs). Intrinsically different galaxy populations exist at high redshifts, which 



cannot be accounted for by these extrapolations. Particularly, Franceschini et al. (2008) 



use observed LFs in the ncar-IR from the local universe to z = 1.4 for describing the ellip- 
tical and spiral populations, and only local for describing irregular/starbursting galaxies. 
They distinguish between these galaxy morphologies using images from different instru- 
ments. Different local LFs and data sets in the IR are used to constrain the mid and 
far-IR background. Their modelling is complex and not reproducible. Despite these par- 
ticular problems, this methodology is based upon LFs, quantity directly observed and 
well understood unlike type (iii) models based on parameterizations of the history of the 
SFR density of the universe, quantity with large uncertainties and biases. 

One important application of the EBL for 7-ray astronomy is to recover the unatten- 
uated spectra of extragalactic sources. This will not be discussed in this proceeding but 
we refer to the interested reader to Dommguez et al. 2011a, b for a discussion about this. 



2. Methodology 



Our model is based on the rest-frame i^-band galaxy LF found in Cirasuolo et al. 



(2010) and on multiwavelength galaxy data from the All- wavelength Extended Groth 
Strip International Survey ( AEGI^ Davis et al. 2 007) of about 6000 galaxies in the 
redshift range of 0.2-1. The Cirasuolo et al. (2010) LF is used to count galaxies (and 



f http://aegis.ucolick.org/ 



EBL from AEGIS galaxy SED-type fractions 
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therefore to normalize the total EBL spectral intensity) at each redshift. The LF as well 
as our galaxy sample is divided into three magnitude bins according to the absolute rest- 
frame K-hand magnitude, i.e., faint, middle and bright. Within every magnitude bin, an 
SED type is statistically attached to each galaxy in the LF assuming SED-type fractions 
that are a function of redshift within those magnitude bins. This is estimated by fitting 
our AEGIS galaxy sample to the 25 galaxy-SED templates from the SWIRE||] library 
(PoUetta et al. 2007). Then, luminosity densities are calculated from these magnitude 
bins from every galaxy population at all wavelengths, and finally all the light at all 
redshifts is added up to get the overall EBL spectrum. 



3. Results and conclusions 

It is shown in Fig. [T]the local EBL, with its uncertainties, compared with direct and 
indirect observational data, and other EBL models. Other quantities such as the EBL 
evolutioijtjare discussed in Dommguez et al. (2011a) Fig.[l]suggests that the EBL coming 
from galaxies is already well constrained in the region from the UV up to the mid-IR but 
not in the far-IR. The EBL measurements free of zodiacal light in two optical bands by 



Matsuoka et al. (2011) agree with our EBL estimations. Furthermore, galaxy counts from 



very deep surveys taken with very sensitive instruments (Madau & Pozzetti 2000; Fazio 
et al. 2004; Keenan et al. 2010) should be considered as a good estimation of the true 
EBL from galaxies. On the other hand, different fully independent modelings based on 
difi'erent methodologies and galaxy data sets such as Franceschini et al. (2008)[ Gilmore 



et al. (2011) Dommguez et al. (2011a) agree in the specific intensity level of the EBL. 



In particular, galaxy count data are in excellent agreement with our EBL estimations. 
From these results, a relevant contribution from Pop III stars to the local EBL seems 
unlikely. 

Summarizing, the best available data sets are used in our modeling (Cirasuolo et al. 
2010's LF and the AEGIS galaxy catalogue) observed over a wide redshift range. This 
model has the following main advantages over other existing EBL models: transpar- 
ent methodology, reproducibility, and -very important- utilizing direct galaxy data. The 
galaxy evolution is directly observed in the rest-frame K band up to z = 4. Observed 
galaxies up to z = 1 from the UV up to 24 /xm with SEDs of 25 different types (from 
quiescent to rapidly star-forming galaxies and including AGN galaxies) are taken into 
account in the same observational framework. A study of the uncertainties to the model 



directly from the data (such as uncertainties in the Schechter parameters of the Cirasuolo 



et al. (2010) LF and the errors in the photometric catalogue) is done. 



It is concluded that the EBL from galaxies seems already well constrained from UV to 
mid-IR wavelengths, even though uncertainties are still large in the far-IR. Furthermore, 
discoveries of 7-ray from distant blazars (e.g., Aleksic et al. 2011a, b,c) support the EBL 
specific intensity level derived from galaxy count and recent EBL models such as Gilmore] 



et al. (2011) Franceschini et al. (2008)^ .Dommguez et al. (2011a) We highlight that the 



EBL specific intensity calculated with our method is matching the lower limits from 
galaxy counts, which implies the highest transparency of the universe to 7-ray allowed 
by standard physics (see Dommguez et al. 2011c). This predicts a promising future for 
the new generation of imaging atmospheric Cherenkov telescopes, namely CTA. 



J http: / / www.iasf-milano.inaf.it /^polletta/templa tes/swire_templates.htm l 
I EBL specific intensities are publicly available at http://side.iaa.es/EBL 
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Figure 1. The solid-black line is the extragalactic background light calculated from our method- 
ology. Direct data, data from galaxy count, upper limits from 7-ray astronomy and other recent 
EBL modelings are shown as well (see Dommguez et al. 2011a for details). 
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